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HOLOMORPHIC TRANSFORMS WITH APPLICATION TO 

AFFINE PROCESSES 

By Denis Belomestny* , Jorg Kampen"!", and John Schoenmakers^ 

Weierstrass Institute 

In a rather general setting of Ito-Levy processes we study a class 
of transforms (Fourier for example) of the state variable of a process 
which are holomorphic in some disc around time zero in the complex 
plane. We show that such transforms are related to a system of ana- 
lytic vectors for the generator of the process, and we state conditions 
which allow for holomorphic extension of these transforms into a strip 
which contains the positive real axis. Based on these extensions we 
develop a functional series expansion of these transforms in terms 
of the constituents of the generator. As application, we show that 
for multidimensional affine Ito-Levy processes with state dependent 
jump part the Fourier transform is holomorphic in a time strip under 
some stationarity conditions, and give log-afnne series representations 
for the transform. 

1. Introduction. Transforms are an important tool in the theory of 
(ordinary and partial) differential equations and in stochastic analysis. In 
probability theory the Fourier transform of a random variable, which repre- 
sents the characteristic function of the corresponding distribution, is widely 
used. Fourier (and Laplace) transforms have become increasingly popular 
in mathematical finance as well. On the one hand, for example, via a com- 
plex Laplace transform and the convolution theorem one may derive pricing 
formulas for one dimensional European options (e.g., see [7]). On the other 
hand, Laplace and Fourier transforms are known in closed form for many 
classes of processes. A famous example is the so-called Levy-Khintchine for- 
mula which provides an explicit expression for the characteristic function 
of a Levy process. More recent financial literature goes well beyond Levy 
processes and attempts to establish explicit or semi-explicit formulas for 
derivatives where underlyings are modelled, for instance, by affine processes 
(see among others [4], [6] and [5]) or more general Ito-Levy processes (e.g. 
see [3]). Theoretical analysis of affine processes is done in the seminal paper 
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by Duffie, Filipovic and Schachermayer [5] and has led to a unique charac- 
terization of affine processes. In particular, it is shown that the problem of 
determining the (conditional) Fourier transform of an affine process X s cor- 
responds to the problem of solving a system of generalized Riccati differential 
equations in the time variable s (see [5]). Although closed form solutions of 
this system can be found in important cases, there is no generic approach 
to solve such a system in the general multi-dimensional case. In this article 
we establish some kind of functional series representation for the Fourier 
transform, hence the characteristic function of the process under consider- 
ation and, in principle, for more general transforms. The most natural one 
is a Taylor expansion in time s around sq = 0. This approach leads to ex- 
pansions where the coefficients can be recursively computed without solving 
linear or non-linear differential equations as, for instance, in the case of the so 
called WKB expansions for transition densities (e.g. see [9]). Unfortunately, 
it turns out that in many cases the resulting power series converges only in 
s = sq = 0. This problem corresponds to a difficulty which is well known in 
semi-group theory and in the theory of parabolic differential equations: small 
time expansions for the solutions of parabolic equations are usually possible 
in a neighborhood of some sq > 0, while an expansion around sq = may 
have zero convergence radius. In this paper we prove that for a generator 
with affine coefficients the Fourier transform extends holomorphically into a 
disc around so = and a strip containing the positive real axes, under some 
mild regularity conditions. Then, for multi-dimensional affine processes we 
obtain convergent expansions for the Fourier transform and its logarithm 
on the whole time line. Hence, we have (affine) series representations for 
the exponent of the characteristic function of a general multi-dimensional 
affine process. More generally, we develop a framework based on a concept 
of analytic vectors which allows for functional series expansions for a class 
of holomorphic transforms which covers the standard Fourier transform. 

The outline of the paper is as follows. The basic setup is described in 
Section 2. In Section 3 we introduce the notion of analytic vectors associ- 
ated with a given generator and study functional series expansions for the 
corresponding transform. Section 4 is devoted to the Cauchy problem for 
affine generators. In Section 5 we derive series representations for the log- 
arithm of the Fourier transform corresponding to a generator with affine 
coefficients. Section 6 gives an explicit representation for these expansions 
in a one-dimensional case. Finally, Section 7 contains results for affine Ito- 
Levy processes which mainly follow from previous sections. More technical 
proofs are given in the Appendix. 



imsart-aap ver. 2007/12/10 file: bks_annal_prob.tex date: July 8, 2008 



HOLOMORPHIC TRANSFORMS 



3 



2. Basic setup. Let (il,^ 7 , (J-t)t>o,P) be a standard filtered probabil- 
ity space where the filtration (J-t) satisfies 'the usual conditions'. On this 
space we consider for each x G W 1 a compensated Poisson random mea- 
sure N(x,dt,dz,u}) = N(x,dt,dz,uj) — v(x,dz)dt on K + x R n , where N is 
a Poisson random measure with (deterministic) intensity kernel of the form 
v(x,dz)dt = E N(x,dt,dz) satisfying v(x,B) < oo for any B G B(W l ) such 
that ^ B (closure of B). Hence TV is determined by 

P[N(x,(0,t],B) =k} = exp(-fa(x, J B)) t - [ X ' B \ k = 0,1,2,... 

A;! 



In particular, for any B G £>(IR) with ^ B and x G M n , the process 
M^' x := N(x,(0,t],B) is a (true) martingale. Further, by assumption, the 
kernel v satisfies, 

v(x,{0})=0, [ (\z\ 2 A \z\)v(x, dz) < oo, i£l n . 



Let us assume W(t) to be a standard Brownian motion in M m living on 
our basic probability space, and consider the Ito-Levy SDE: 

(2.1) dX t = b(X t )dt + a{X t )dW{t)+ [ zN{X t -,dt,dz), X = x, 

for deterministic functions b : R n — > W n , a : K n — > M n x R' m , which satisfy 
sufficient regularity and/or mutual consistency conditions such that (2.1) 
has a unique strong solution t — > Xt, called an Ito-Levy process, which can 
be regarded as a strong Markov process (e.g., see [11], [3]). 

As a well-known fact, the above process X can be connected to some kind 
of evolution equation in a natural way. In this context we consider a 'pseudo 
generator' 

(2.2) A* : V(A^) c C (2) cC — > C, 

where C := C(M n ) is the space of continuous functions / : M n — > C, equipped 
with the topology of uniform convergence on compacta, and is the space 
of functions / G C which are two times continuously differentiable. Further, 
/ G V{A^) iff / G and 



(2.3) A* f(x) := - £ ^(*)^+£^)^ 

w(x, cfe), with a := <7<7 T , 



2 .•4-' 1 dxidxi r-f dxi 

1,3=1 ■> 1=1 



+ 



df 

f{x + z) - f(x) - — ■ z 
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exists and is such that A*f E C. In this respect we assume that the building 
blocks b(x), a(x), v(x,dz) of SDE (2.1) have bounded derivatives of any 
order with respect to x. Clearly, T>(A^) is dense in C and it can be shown 
that the operator A$ thus defined is closable. The closure of A^ is denoted 
by 

(2.4) A : V(A) cC — ► C. 

As such A can be seen as a relaxation of the notion of a generator of a 
strongly continuous Feller-Dynkin semigroup associated with the process 
X, for which the Hille-Yosida theorem applies. This semigroup is usually 
defined on the Banach space Cq(W 1 ), i.e. the set of continuous functions 
on W 1 which vanish at infinity, equipped with the supremum norm, and its 
generator coincides literally with A^ on a dense subdomain of Co(M n ). By 
slight abuse of terminology however, we will also refer to A as 'generator' 
when A is considered in connection with the process X given via (2.1). Let 
now $ '■= {fui u £ 1} C C be a dense subset of bounded continuous functions 
f u : W 1 — > C which have bounded derivatives of any order. With respect to 
the (closed) generator (2.4) we consider for each f u € 5 the (generalized) 
Cauchy problem 

dp 

(2.5) —(s,x,u) = Ap(s,x,u), 

p(0,x,u) = f u (x), s>0, xeJcK", 

where X is some open (maximal) domain, and assume that problem (2.5) 
has a unique solution. For instance, see [1] for mild conditions which guar- 
antee existence of unique global solutions of integro-differential evolution 
problems. In particular, if some ellipticity condition is satisfied we may have 
X = W 1 . For mixed type generators, such as affine generators, existence, 
uniqueness, and the maximal domain X has to be considered case by case. 
For example, if b(x), (cra T )(x), v(x, dz) are affine in x and satisfy the admis- 
sibility conditions in [5], existence and uniqueness are ensured in a domain 
of the form X = M 1 x c K n (for details see [5]). We underline, however, 
that in this article the main focus is on functional series representations for 
the solution of (2.5), and we therefore merely assume that sufficient regu- 
larity conditions for the coefficients in (2.3) (hence (2.1)) are fulfilled. 

Remark 2.1. In our analysis we often consider the pseudo generator 
(2.3) and it's closure (2.4) on C := C(X), for an open domain X C W 1 , rather 
than C(M n ). For notational convenience (while slightly abusing notation) we 
will denote these respective operators with A^ and A also. 
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If A is the generator of the process (2.1), the solution p(s,x,u) has the 
probabilistic representation 



0.2N 



p(s,x,u) = E f u (X_ 



where X 0,x is the unique strong solution of (2.1) with Xq' x = x a.s. We refer 
to p(s, x, u) as generalized transform of the process X® ,x associated with 
As a canonical example we may consider 



(2.6) 



f u (x) := e 



in which case (2.5) yields the characteristic function p(s, x, u) = E[e m Xs ' ]. 

By using multi-index notation, the integral term in (2.5) may be formally 
expanded as 

df 1 

f(x + z) - f{x) - — ■ z v(x, dz) = 

If 1 

X] ~^] dxa f( x } / z0lv { x 'dz) =: X —m a (x)d xa f. 



\a\>2 



\a\>2 



Hence, we may write formally the generator as an infinite order differential 
operator 



(2.7) 



A = a a (x)d xa 

\a\>0 



with obvious definitions of the coefficients a a (x) for a > 0. 

3. Analytic vectors and transforms. First we introduce the notion 
of a set of analytic vectors associated with an operator A. 

Definition 3.1. 5 = {/ u , u € /} is a set of analytic vectors for an 
operator A in an open region X, if 

(i) A k f u exists for any u £ I and k G N, 

(ii) for every u E I there exists R u > such that for all x G X, 



E |r ^(*) 



k=0 



< oo, < s < R u , 



where the convergence is uniform in x over any compact subset of X. 
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Thus, if ^ is a set of analytic vectors in the sense of Definition 3.1 then 
for all x G X the map 

00 s k 

(3.1) s P s f u (x) := ]T T] A k f u {x), \s\ < R u 

k=o K - 

is holomorphic in the complex disc Dq := {s € C : \s\ < R u } and the series 
converges uniformly in x over any compact subset of X. In fact, P s f u (x) 
coincides with p(s, x, u) for < s < R u . 

Proposition 3.2. If $ is a set of analytic vectors, the map (s,x) — > 
Psfu{%) defined in (3.1) satisfies (2.5) for all s, \s\ < R u and x £ X. In 
particular we have P s f u { x ) = p( s i x ,u), < s < R u . 

Proof. Obviously, Pof u (x) = f u (x) for x € X. Set (see Remark 2.1) 

N k 

P s {N) fu(x) := £ yA k f u (x), 

k=0 

then both Pg N ^ f u (x) and 

N k 

AP^Ux) := £ ^U fc+1 / U (*) 
fc=o K - 

converge uniformly for any a; in a compact subset of X and for any s satisfying 
|s| < R u — e with arbitrary small e. Hence, since A is closed, 

00 k pi 00 k o 

,4P S / U (x) = £ ^U fe+1 / U (z) = ^ £ lr-A*/»(x) = «-P./«(x). 



□ 



In order to study generalized transforms associated with a set of analytical 
vectors 5 hi domains containing the non-negative real axis we introduce for 
n > the sequence 



(3.2) qi>\x,u) := ^ll^A + r^U 



k-l 



k\ 

r=0 

k 



X) 



(3.3) =:±-S2c k , r r,- r A r f u (x), x e X, u€l, fc = 0,1,2, ... 



A;! 

r=0 



imsart-aap ver. 2007/12/10 file: bks_annal_prob.tex date: July 8, 2008 



HOLOMORPHIC TRANSFORMS 



7 



In (3.3) the coefficients C)., r , < r < k, are determined by the identity 

k-1 k 

11 [z + r) = z(z + 1) • . . . • (z + k - 1) = c k,rz r , 

r=0 r=0 

and are usually called unsigned Stirling numbers of the first kind. These 
numbers satisfy co,o = 1 and 

c fe,o = 0, Ck,k = 1> 

(3.4) c k+1 j = kc k ,j + c k)j -i, l<j<k, 

if k > 1. Obviously, the following recursion is equivalent to (3.2), 

(3.5) (k + ^q^fau) =r ] - 1 Aq ( k r > ) (x,u)+kq ( j? ) {x,u), k>0, x€l 

The next theorem provides a functional series representation for the solution 
of (2.5) for all s > 0, under certain conditions. 

Theorem 3.3. Let $ be a set of analytic vectors in the sense of Defini- 
tion 3.1, u € / be fixed, and the sequence (<zl ) be defined as in (3.3). Let p 
be the solution of the Cauchy problem (2.5). Then the following statements 
are equivalent: 

(i) There exists a constant R u > such that for each x E X, the map 
s — ► p(s,x,u) has a holomorphic extension to the domain 

Gr u ■■= {z : \z\ < R u } U {z : Rez > A | Imz\ < Ru} , 

see Figure 1. 

(ii) There exists an r] u > such that for each x € X the following series 
representation holds: 

oo 

p(s,x,u) = Y / q < ^ u) (x,u)(l-e- VuS ) k , < s < oo. 
fc=0 

Moreover, the series converges uniformly for (x, s) running through any com- 
pact subset of X x {s £ 1 : < s < R u } . 

(Hi) The solution p of the Cauchy problem (2.5) is holomorphically ex- 
tendable to [0,oo), there exists r] u > such that 



(3.6) limk^ {I q^\x,u) < 1, x G X, 
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Fig 1. Domain Gr u on the complex plane 



and, there exists e u , < e u < 1, such that the series 



(3.7) 



fc=0 



u)w 



converges uniformly for (x, w) running through any compact subset of X x 
{w G C : \w\ < 1 — e u } . 



Proof. See Appendix. 



□ 



Remark 3.4. The implication (iii)' =>■ (i), where statement (iii)' con- 
sists of (3.6), and (3.7) with e u = holds as well. That is, loosely speaking, 
if in (iii) series (3.7) converges uniformly on all compact subsets of X x 
{w £ C : \w\ < 1} , the holomorphy assumption on p can be dropped. 

Remark 3.5. In order to use the representation in (ii) one has to choose 
rj u . In fact, rj u can be related to R u via r/ u = tt/R u and hence increases with 
decreasing R u . 

It is important to note Theorem 3.3 concerns the solution of the Cauchy 
problem (2.5) connected with a general operator A. In particular, all criteria 
in this theorem are of pure analytic nature and via (3.2), respectively (3.3), 
exclusively formulated in terms of the A k f u (x), i.e. coefficients in Defini- 
tion 3.1. In the case where A is the generator of a Feller Dynkin process one 
can formulate a sufficient probabilistic criterion for Theorem 3.3-(i): 
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Proposition 3.6. Let$ be a set of analytic vectors in the sense of Def- 
inition 3.1 and let the Markov process {Xt} be associated with the generator 
A. If in addition, for every u € / there exists a radius R u such that for any 
t > 



00 „fc 

( 3 - 8 ) E it Ak nfu(x?n 



< oo, < s < R u , 



uniformly in x over any compact subset of X, then Theorem 3.3-{i) holds. 

The statement is a direct consequence of the following "quasi" semi-group 
property of the transition operator P t . 

Proposition 3.7. Let $ be a set of analytic vectors satisfying (3.8). 
Then, for all x G X and all t > 0, the generalized transform p(t + s, x, u) 
can be represented as 

(3.9) p(t + s,x,u) = J2 T7A k E[f u (X?' x )], 0<s<R u , 

k=o fc - 

where the series converges uniformly in x over any compact subset of X. 

Proof. Denote the right-hand-side of (3.9) by p(t, s, x,u). Obviously, 
p(t,0,x,u) =E[f u (X?' x )}. Set 

N k 

p^(t, S ,x,u) :=Y. S -nA k nUX^)\, 
k=o K - 

then both (t, s, x, u) and 

N k 

A^ N \t,s,x,u) = £ S -A k ^nUX Q t ' x )] 
k=o K - 

converge for N — > oo uniformly over any compact subset of X, and s with 
|s| < R u — e, for an arbitrary small e > 0. Hence, for |s| < R u — e, we have 

a N-X k 

gpW(t, s ,x,u) = £ S -A k + 1 M[f u (X^)]=Ap< N - 1 \t, S ,x,u), 

k=0 

p(t, 0, x, u) = p(t, x, u) 

and thus, by closeness of the operator A and uniqueness of the Cauchy 
problem (2.5)-(2.7), we have p(t, s, x, u) = p(t + s,x,u). □ 
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The following proposition provides a situation in a semigroup context 
where a much stronger version of the condition (i) in Theorem 3.3 applies. 
It also sheds light on the connection between semi-group theory and holo- 
morphic properties of generalized transforms. 

Proposition 3.8. Let Co(M. n ) be the Banach space of continuous func- 
tions f : W 1 — > C which vanish at infinity, equipped with supremum norm: 
ll/H : = sup^n \f(x)\. Let A : 33(A) C C (R n ) -► C (M n ) be the generator 
of the Feller-Dynkin semi-group (P s ) s >q associated with the process X, i.e. 
P s f{x) = E[f(X°> x )], f G Co(M n ). Suppose that the family £ is such that 
f u G 33(A k ) for each u G I and all integer k > 0, and that for each u £ I, 



E ¥ 



< oo, < s < 



T/ien /or eoc/i u £ L, 

(3.10) P s / U = ]T < s < R u , 



oo gk 



n kl 

k=0 



with convergence in Co(M n ). Thus, the map s — > P s f u for < s < R u extends 
via (3.10) to the complex disc Dq := {s € C : \s\ < R u } . In particular, for 
each x G M. n the map s — > P s f u (x) is holomorphic in Dq. Moreover, for each 
t > 0, we may extend the map s — > Pt+ S fu, < s < R u to the disc Dq via, 

00 „fc 00 „k 

(3.11) P s+t f u = P t P s f u = -H P t Ak fu = E ~u Akp tfu, s G D . 

k=o fc - fc=o K - 

Proof. See Appendix. □ 

Thus, under the conditions of Proposition 3.8, £ = {f u , u G /} is a set 
of analytic vectors for the generator A in the sense of Definition 3.1 with 
X = K n . Moreover, due to Proposition 3.8 the map 

oo k °° k 

s - P s+t f u (x) = £ S -P t A k f u {x) = £ i^Et/upr?*)], 
fe=o K - fc=o K - 

is holomorphic in Dq for each x G K n and hence Theorem 3.3-(i) is fulfilled. 

In this paper we do not stick to the semigroup framework because we 
want to avoid the narrow corset conditions of Proposition 3.8. We also want 
to consider operators A with unbounded (for instance, affine) coefficients 
and sets 5 of functions that do not vanish at infinity (for example, (2.6)). 
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Such situations may lead to the violation of condition f u € D{A k ), k € N 
in the sense of Proposition 3.8. In particular, in the next Sections 4-5 we 
will focus on general operators A with afhne coefficients and in Section 7 on 
affine processes related to affine generators satisfying a kind of admissibility 
conditions. 

4. Affine generators. Let us now consider generators of the form (2.3) 
with affine coefficients. It is important to note that in this section A may or 
may not be a generator of some Feller-Dynkin process. The next theorem 
and its corollaries say that the Fourier basis ^ = {fu( x ) '■= e m x , u G M n } 
is a set of analytical vectors for an operator A of the form (2.7), where 
the coefficients a a (x) are affine functions of x and satisfy certain growth 
conditions for \a\ — » oo. Moreover, an explicit estimate for the radius of 
convergence is given. 

Theorem 4.1. Let A be a generator of the form (2.7) with affine coef- 
ficients a a (x), i.e. for all multi-indexes a, 

(4.1) o. a {x) =: c a + x T d a , x 6 X, 

where X C W l is an open region, c a is a scalar constant, and d a £ M. n is a 
constant vector. Assume that the series 

(4.2) ]TM*)(iur 

H>o 

converges absolutely for all u € M. n . Then, for every u S M n and x E X it 
holds 

(4.3) \A r f u (x)\<(r + l)\2 nr (l + \\x\\ye r (\\u\\) 

with \\x\\ = max \xi\, 
i=l,...,n 

(4.4) 9(v) :=J22 k (l + v) k V a k , v£R + 

k>l 

\d xf >a a (x)\ 
sup max — — . 

xeX \a\=k,\p\<i l + \\x\\ 

The proof of Theorem 4.1 is given in the Appendix. 

Corollary 4.2. If in Theorem 4-1 the region X is bounded, the Fourier 
basis J constitutes a set of analytic vectors for the affine operator A in X. 



and 

n ■= 
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Corollary 4.3. If in Theorem 4-1 there exists for any ? > a constant 
M ( which may depend on ^ > 0) such that 

V a k < Mq k /k\ k > 1, 

then 

9(v) < Afexp(2?(l +«)) . 

Corollary 4.4. // in Theorem 4-1 it holds that a a (x) = for \a\ > 2 
( generator of diffusion type ) then 

0(v) < C(l + v) 2 , C>0. 

Remark 4.5. The requirement that (4.2) converges for all u imposes 
restrictions only on the tails of the measure v and so does not exclude 
infinite activity processes. 

For an affine operator A the sequence (3.2) can be explicitly constructed 
via the next proposition, which is proved in the Appendix. 

Proposition 4.6. Let A be an affine generator as in Theorem 4-1 and 
define 

bp(x,u) := d u0 = 2^ a a+f3 (x) - (m) 

Ju{X) a > Ct. 
K, |k| = i 

with do := 0. We set A r f u (x) =: g r (x,u)f u (x) and, for fixed rj > 0, qi v \x,u) =: 
h r (x, u)f u (x), where both g r and h r are polynomials in x of degree r. It holds 

(4.5) g r (x,u)=: ^ gr^^x 1 , h r (x,u) =: ^ K^^ujx 1 , 

\j\<r | -y | < r- 
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where g r and h r satisfy go = <?o,o = ^0 = ^0,0 = 1) and for r > 0, respectively, 
(4.6) g r+hl = \ + p ]9r,y+l3bp 



\@\<r- |t| 



+ J2 2 ( 7 ^ + /3 j5r, 7 - K +/3^, K ' alld ' 



|k|=1,«<7 |/9|<r+l-|7| 



\/3\<r-\~f\ \ 1 ) 



| K |=1, K <7 |/3|<r+l-|7| \ P / 

where \^\ < r + 1, and empty sums are defined to be zero. 

Remark 4.7. Depending on the open set X we may consider instead of 
(4.5) for an xq € X expansions in x — x$ rather than in x. For simplicity 
we henceforth assume {0} € X which, if necessary, may be realized by a 
translation of the state space. 

A natural question is whether affine generators are the only ones for which 
the Fourier basis constitutes a set of analytic vectors. For this paper we leave 
this issue as an open problem but the following proposition shows that at 
any case the set of such generators is rather "thin". 

Let us put X = [—it, n] and 

(4.7) A= l aW ^ +6(x) ^. 

Proposition 4.8. The set of coefficients (a(x),b(x)) such that for an 
arbitrary M > 

\\A N f u \\ L 2 (x) > M N Nl , iV^oo, 
is dense in L 2 (X) x L 2 (X). 

Proof. Without loss of generality let us assume that b(x) = and u > 0. 
The general case can be considered along the same ideas and is only formally 
more complicated. Any a £ L 2 {X) may be approximated (in L 2 -sense) by a 
finite Fourier series 



(4.8) a(x)^J2 a i eilX - 

1=1 
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Thus, for given e > we can find natural n and amplitudes ai [a n ^ 0) such 
that 



(4.9) 



i=i 



tlx 



< £. 



The corresponding approximative operator is given by 
(4.10) 



l=i l=i ox 



Using the fact that for any s\, . . . , s& € N, 

fc-i 



A Sl -..A Sk e- = (-l) k a sl .--a Sk ]l « + 5> ) e^^, 

1=0 \ j=0 



and setting 



F k := 



1 



e- ikx f-\x)A N U{x)dx, ken, 



2tt J- w 

we have Fk = for k > nN, and for N — ► oo 
N-l 



F nN = (-lfa» l[(u + nlf ~ (-l«n 2A? ((iV - l)!) 2 (iV - I) 2 "/™. 



1=0 



Further, by Parseval's identity it holds 



A N fu\\L 2 (X) 



nN 



.1/2 



fc=0 



and then we are done. 



□ 



Obviously, Proposition 4.8 may be formulated with respect to any compact 
interval. 

5. Log-affine representations for the afflne Cauchy problem. Let 

us consider the Cauchy problem (2.5) for affine generators A of the form 
(2.7), under the assumption (4.2). As in (4.1) we set a(x) = c a + x T d a . The 
ansatz 



(5.1) 



p(s, x, u) = exp (c(s, u) + x T D(s, uu , 
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for scalar C(s, u) and vector valued D(s, u), where C(0, u) = and D(0, u) = 
iu, for the Cauchy problem (2.5) yields, 

d s C + x T d s D= Y a a (x)D a = Y c a D a + Y x T d a D a , 

\a\>0 |ct|>0 |q|>0 

and so 

8 S C = Y c «D a , d s D = Y, d aD a , 

\a\>0 \a\>0 

C(0,u)=0, D(0,u)=xu. 

We thus have a system of ordinary differential equations (ODEs), which 
reads component- wise 

d s C = Y c * Da > d s Dj = Y d^D a , J = 

q|>0 |a|>0 

(5.2) C(0,n)=0, Dj(0,u)=Wj. 
By assumption (4.2), the series 

(5.3) y c «y a > Y. d( i ] y a i j = i,-.«> 

|q|>0 |a|>0 

are absolutely convergent for all y G M n , and thus define terms- wise dif- 
ferentiable C^°°)(R n ) functions. In particular, they are locally Lipschitz 
and so according to standard ODE theory the system (5.2) has for fixed 
u £ W 1 a unique solution (C(s,u), D(s,u)) for < s < < oo, where 
(s, C(s, u),D(s, u)) leaves any compact subset of R x M x W 1 , when s | sjf. 

Remark 5.1. By a general theorem from analysis (e.g., see [2]), it follows 
that the solution of (5.2) extends component- wise holomorphically in s into 
a disc around s = 0, due to the analyticity of (5.3). This implies that (5.1) 
is holomorphic in s. So, besides Theorem 4.1, also along this line one may 
show that (5.1) can be represented as a power series of the form (3.1). I.e., in 
particular, the Fourier basis (2.6) constitutes a set of analytic vectors for the 
affme operator A. However, the direct approach in the proof of Theorem 4.1 
(see Appendix) leads to an explicit estimate (4.3) and allows for investigating 
possible extensions of p(s, x, u) into a strip containing the real axis in the 
complex plane (see Theorem 5.4). Moreover, it also suggests the line to follow 
in cases where A is not affine and/or the function base is not of the form 
(2.6). 
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Let us suppose that for fixed u £ M. n the statements of Theorem 3.3 hold. 
Then we obtain for < s < < oo, x € X, 



(5.4) 



p(s, x, u) = exp (c(s, u) + x T D(s, u) 



k=0 



Since q$\x, u) = f u (x) = exp (iu T x^ ^ we have, taking into account the 
boundary conditions for C and D, at least for small enough e > 0, 



(5.5) C(s,u) + x T D(s,u) = J2p^ u) (x,u)(l -e~^ s ) k , < s < e, 

k=0 

where by a standard lemma on the power series expansion of the logarithm 
of a power series, 



(5.6) 



p { o u \x,u) = lnq ( Q lu) (x,u) = iu ' x, 



,07"), 



fu{x) 



k > 1. 



Thus by (5.5), the p£ u ^ are necessarily affine in xl 

(v ) 

Remark 5.2. It is possible to prove directly that the functions pY 
defined above are affine in x using Proposition 4.6 via a (rather laborious) 
induction procedure, so without using a local solution of (5.2). 

Theorem 5.3. Suppose that for fixed u £ I the statement of Theo- 
rem 3.3- (i) holds for an open region X and, in addition, for s € Gr u and x € 
X it holds p(s,x,u) + 0. Then, with p^ u \x,u) =: p% u '°\u) + x T p^ (u) 
determined by (5.6), we have 



p{s, x, u) = exp 



oo 



lk=0 



< s < oo. 



Proof. Let u G I and x € X be fixed. Since Gr u is simply connected 
and s — > p(s, x, u) is holomorphic and non-zero in Gr u , there exists a branch 
s — ► L(s,x,u) of the logarithm such that p(s,x,u) = exp(L(s, x, u)) for 
all s € G/j u . Along the same line as in Theorem 3.3 we then argue that 
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there exists an r\ u > such that w — > L($> Vu (w),x,u) (see the proof of 
Theorem 3.3) is holomorphic in the unit disc {w :| w |< 1}, hence, there 
exists pk(x,u) such that 

L(<k Vu (w), x, u) = ^ pk(x, u)w k , < | to | < 1, and so 

fc>0 

L(s,x,u) = ^2 Pk(x,u)(l - e~ VuS ) k , < s < oo. 

fc>0 

Since the later expansion must coincide with (5.5) for small s, it follows that 
necessarily pk(x,u) = p^ u \x,u) and the theorem is proved. □ 

Let us now pass to another interesting log-afnne representation for the 
characteristic function. Prom (4.5) and (4.6) we derive formally 

5>W(s,u)(l - e-"T = e iuTx Yl E h r , 7 (u)x^(l - e^Tl,^ 

r>0 r>0 7>0 

7>0 r>0 

Suppose that the requirements of Theorem 5.3 hold. Then, using Theo- 
rem 4.1, it is easy to show that for small enough e > 0, 

E H X Hoo E l /l |7|+r,7( n )ll u 'l' 7 ' +r < °°) if H < e ' IHIoc < £ 
7>0 r>0 

(see Remark 4.7). Thus, for |s| and || 

x|joo sm&U enough we obtain 



lnp(s,x,u) = in T x + ln ^x 7 ^/i| 7 | +rj7 (u)(l -e-^)! 7 ^' 

\7>0 r>0 

= C(s, u) + x T D(s, u), 
with (in multi-index notation) 



(5.7) C(s, u) = In \J2 M^X 1 - e~^ s 

\r>0 

D (s,u)=m + - ~ - , , n — — , |«| = 1. 



Er>0^,o(«)(l-e-^) f 
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However, the left- and right-hand-sides of (5.7) are holomorphic for all s € 
Gr u and we so arrive at the representation 



(5.8) p(s,x,u) = exp 



In ( X>r,oM(l -e-<* a y] +iu T x 

r>0 



T Er>l Mu) (1 - e-"» l 

X 



, s e G Ru) x g X, 



Er>cA,o(«)(l-e-^ s ) r 
with 

&r(u) := Ki(")] i=lv .. jn , 
where for 1 < i < n, the multi-index (^tj)j=i,---j« * s identified with i. 

Particularly due to the explicit estimate (4.3) for affine generators in 
Theorem 4.1, we may proof the next theorem (which is a non-probabilistic 
version of Proposition 3.6 in the situation where A is affine). 

Theorem 5.4. Let X be a bounded domain. Assume that the system (5.2) 
is non- exploding, i.e. s^ 3 = oo, and that for any fixed u £ M. n the solution 
D(s,u) remains bounded as s — > oo. Then, there exists a radius R u > suc/i 
i/iai for any t > f/ie map s — > p(i + s, x, u), < s < R u has a holomorphic 
extension to the disc {s € C : |s| < R u }- Moreover, it holds 

°° s k 

(5.9) p{t + s,x,u) = ^—A k p{t,-,u){x), \s\ <Ru, x e X. 

Remark 5.5. The maximal extension radius R u satisfies 

(5.10) R u > — T-777— ; .... inf - — 77, 

2 n e(\\D*(u)\\) xexi + \\x\y 

where function 6 is defined in (4.4) and ||Z)*(«)|| = sup s>0 \\D(s, u)\\. 

Proof. Denote the right-hand-side of (5.9) by p(t, s, x,u). Obviously, 
p(t,0,x,u) = p(t, 0, x, u). Let us define 

N k 

p( N \t,s,x,u) : -=^2^A k p(t,;u) (x). 

Since p(t,x,u) = exp(C(t,u) + x T D(t,u)), Theorem 4.1 implies that the 
series 

00 

(5.11) Ett^V^H*) 

fc=o K - 
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is absolutely and uniformly convergent on any compact subset of X x {s € 
C : \s\ < R u } if R u satisfies (5.10). So, both p^ N \t, s, x, u) and 

ApW(t,s,x,u) = J2TrA k+1 p(t,;u) (X) 

converge for TV — > co uniformly over any compact subset of X and s with 
|s| < R u — e for any e > 0. Hence, for \s\ < R u — e 

ft N-l k 

Q-/ N \t,s,x,u) = £ S ^A k+1 p(t,;u) = Ap( N - l \t,s,x,u), 

p(t,0,x,u) =p(t,x,u), 

and thus, by closeness of the operator A and uniqueness of the Cauchy 
problem (2.5)-(2.7), we have p(t, s, x, u) = p(t + s,x,u). □ 



6. Full expansion of a specially structured one dimensional affine 
system. Let us consider Cauchy problem (2.5) for n = 1 with f u (x) = 
exp(inx), where the jump-kernel in the generator A (see (2.3)) has a special 
affine structure of the form 

v(x, dz) =: (Ao + Xix)n(dz), 

and where the diffusion coefficients have a similar structure, 

b(x) = (Ao + Xix)9, a(x) = (Ao + Aix)t9, 

for some constants Ao,Ai,0,i9 G K, and measure fi. So, in Proposition 4.6 
the a a have the form a/ =: (Ao + X\x)r]i where 

r) := 0, ??! = 9, r] 2 :=-U + j z 2 fi(dz)^j , ru := - J z l fi{dz), I > 2. 

Hence, in Proposition 4.6, the bp in (4.6) have the form 
b r (x,u) = br(u) + xbl(u) 

=: (Ao + Xix) m+r^jf^iiu) 1 =: (A + \ lX )^(u) r > 0, 
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where f)(u) := J2i>o r 1i(^ u ) 1 - It is now possible to show via (4.6) that for 
r > 1, 



P>0,q>0 
0<ni<...<n q , mi,...,m g >0, 
r=p+mmiA \-n q m q 



X 

i=i 



with the following integer recursion procedure: 

Initialization: ttGO = 1, ^ ( ( ° ) limi) ,..., K , mg) =0, p, ? > 1. 
For all nj > 0, mi > 0, with 1 < i < p, g > 1: 
Reduction rule I: If rrij = 0, for some j, 1 < j < q, then 

-GO = -GO 

(ni,mi),-,(n„m,) (m,TOi),-",(rty-i,m^_i),(ry+i,?Ty+i),"- ,(n g ,m 9 ) " 
Reduction rule II: 

» 

q 



7T 

(ni,mi),-,(n,_i,m,-i),(n„m,) 



^,p + nj -l\ (p+ry-1) W P_1) 
^— ^ ^ / (ni,'mi),---,(nj,mj-l),---,(n g ,m ? ) (ni , mi ), ■■■,(«,, m 9 ) " 



In fact, the above recursion procedure follows automatically after substitut- 
ing (6.1) as ansatz into (4.6). 
Note that 

(«*)' f J. 



&(u)=i«e-^ 2 + X;^f-/« l M(d 



i>2 

i)-l + (i#-0'(O))u-^u 2 , 



where ^ is the characteristic function of the measure \x. Hence, f) and all its 

(v) 

derivatives may be computed from (p. Subsequently we obtain the qi for 
the series expansion in Theorem 3.3 by (3.3), i.e. 



q { k\x,u) y^2c k , r 7] r g r 



k 

L 

r=0 
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7. Application to affine processes. Affine processes have become 
very popular in recent years due to their analytical tractability in the context 
of option pricing, and their rather rich dynamics. Many well-known models 
such as Heston and Bates models fall into the class of affine jump diffusions. 
Option pricing in these models is usually done via the Fourier method which 
requires knowledge of the Fourier transform of the process in closed form (see 
e.g. [4]). The functional series representations for affine generators developed 
in this paper, in particular (5.8), can be directly applied to affine processes. 
Let us recall the characterization of a regular affine process as given in [5]. 

Definition 7.1. We call a strong Markov process {Xt} with generator 
A a regular affine process if A is of the form (2.3) and all functions 

a,ij(x), bi(x), v(x,dz) i,j = l, ...,m 

are affine in x (see (4.1)), and satisfy the set of admissability conditions 
spelled out in [5, Definition 2.6]. These conditions guarantee that A is the 
generator of a Feller-Dynkin (strong) Markov process A in a subspace of 
the form R l x R^~ l C R n for some < / < n. 

The next theorem provides a sufficient condition for convergence of the 
series representation in Theorem 3.3-(ii), hence representation (5.8), for reg- 
ular affine processes. 

Theorem 7.2. Let {X s } be a regular affine process which has a non- 
degenerated limiting distribution for s — > oo, and has a generator A which 
satisfies the moment condition (4-2). Then the (conditional) characteristic 
function p(s , x , u) = E,{f u (X® ,x )], with f u (x) = e m x , has a representation 
according to Theorem 3.3- (ii): 

oo 

p(s, x ,u) = J2 Qk u \ x i «) (! - e~ VuS ) k , < s < cxo. 

k=0 

Moreover, the scaling factor r] u may be chosen according to the inequality: 

Vu >ce(L(i + \\u\\ 2 2 )), 

where the (monotonic) function 6 is defined in (4.4), L > is a constant 
independent of x, and C is a constant generally depending on x. 

Proof. Following [5], p(s,x,u) has representation of the form (5.1) for 
< s < oo. The existence of a limiting distribution implies in particular 
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that D(s,u) in (5.1) is bounded for all s > 0. Moreover, as shown in [5, 
Section 7], p(s,x,u) is the characteristic function of some infinitely divisible 
distribution for all s > (hence also in the limit s — > oo). As a consequence 
(see [12]), there exists an M > independent of s such that 

lim IMI2 2 |l°gp( s i x i u )\ < M. 
||m||2— +00 

This implies that ||D(s,n)|| < L(l + for some constant L > not 

depending on x and s > 0. Now we apply Theorem 5.4 and Remark 5.5. □ 

Remark 7.3. The existence of a limiting stationary distribution is a 
sufficient condition for the boundedness of D(s,u). In fact, there are affine 
processes which have no limit distribution but bounded D(s,u) (a trivial 
example is standard Brownian motion). The study of existence of limiting 
(and stationary) distributions for affine processes is currently under active 
research, e.g. see [10] or [8]. 

8. Appendix. 

Proof of Theorem 3.3. (i) =^> (ii) : Let U := {z G C : \z\ < 1} be the 
unit disc in the complex plane. Consider for rj > the map 

: z — ► — Ln(l - z). 
rj 

Obviously, there exists an rj u > such that 

(0- co)c<MW)cG^. 

Moreover, the map $„ u is injective on IA. Thus, (denoting the extension in (i) 
with p as well) for each the function p($ Vu (w), x, u) is holomorphic 

in hi and has a series expansion 

00 

w — >p(<Z> Vu (w),x,u) =: qk{x,u)w k , \w\ < 1, 

fc=0 

and as a consequence, 

00 

(8.1) p(z,x,u) = Y j qk{x,u){l-e-^ z f, z&$ Vu (U), x G X. 

k=0 

Since (8.1) holds in particular for z G (0—, 00), we have in a (possibly small) 
e-disk around z = 0, 

00 00 b 

(8.2) p(z,i, U ) = ^g fc (x lU )(l-e-"»f = ^-i fc /uW, < \z\ < e, 

k=0 k=0 
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due to Proposition 3.2. By taking z = we have 

p(0,x,u) = q (x,u) = f u (x). 
Taking derivatives at z = yields 



dz k 



k oo 



1=0 



k 3 k I /A 

= E^)|,E H) v 

=0 1=0 3=0 \ J J 



Z = 



1 = 7=0 V/ 



2=0 



hence 
(8.3) 



1=0 3=0 



J, 



-iy(-j) k = vu k A k fu(x). 



In Lemma 8.1 it is proved that the solution of system (8.3) is given by 
qfau) = q\ nu \x,u) with q^ defined in (3.3). 

Next we prove the uniform convergence as stated in (ii). A well known 
property of Stirling numbers implies that the series 

\k oo 



E \vu r AT fu{x)\ J2 c k,r-^r = E 



r=l 



k=r 



r=l 



A r f u (x) 



V u r |!og(l - \w\ 



converges uniformly on any compact subset of X x {w : \w\ < 1 — e 
Thus, by a Fubini argument the series 

00 n _ P -vus\k k oo 

E 1 \ ) E c fe ^-M74^) = E <it\*, «)(i - 



fc=0 



r=l 



fc=0 



is also uniformly convergent on any compact subset of X x {s : < s < R u } . 
(ii) => (Hi) : Is obvious, take e u := 1 — e~ nuRu . 

(Hi) ==> (i) Let Tj u and e u be such that (m) holds. We may then define 
(see the proof of (ii)) 



(8.4) 



p(z,x,u) = J^q^ u) (x,u)(l -e~^ z )\ x G X, 
k=0 



which is holomorphic in z G & Vu (U). We first note that p(0,x,u) = f u (x)- 
Next we consider for ^ s <C — 7] u ^ ln£ w , 



p< N Hs,x,u) = Y,qt\x,u)(l- e 

k=0 



■r] u s\k 
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which satisfies 

N 







p^ N \s,x,u) = ]T kq [ ™(x,u)(l - e^f-V""" 
as k=l 

N 



k=l 
N-l 

(8.5) +rju J2(k + l)q^l(x,u)(l-e-^ s ) k 

k=0 

= -r] u Nq^\x,u)(l-e-^ s ) N 
N-l 

(8.6) + J2Aq^\x,u)(l-e-^ s ) k 

k=0 

by some rearranging and using (3.5). Since due to (Hi) the first term in (8.6) 
vanishes for N — > oo, we obtain 

-^-p(s,x,u) = lim -^-p^ N \s, x, u) = lim Ap^ N ~ 1 \s, x, u), 

OS N-*oo OS N-^oo 

together with 

lim p( N \s,x,u) = p(s,x,u). 

N-+oo 

From the uniform convergence as stated in (iii) it follows easily that the two 
series in (8.5), the first term in (8.6), and so also the second term in (8.6) 
converge uniformly in the same sense. Thus, the above limits are uniform 
on compacta accordingly. Since the operator A is closed, we so obtain 

q _ 

—p(s,x,u) = Ap(s,x,u), < s < -rj~ lne«, 

OS 

and by uniqueness of the Cauchy problem associated with the operator A 
we thus have 

N 

p(s,x,u) =p(s,x,u) = J2lk u) ( x ' u )( 1 -e~ VuS ) k , < s < -Ti^lne^. 

k=0 

Because of the assumption that p(s,x,u) is holomorphically extendable in 
each s, < s < oo, we then must have p(s, x, u) = p(s, x, u) for < s < oo. 
Finally, it is not difficult to see that there exists R' u > such that Gr> u C 
<& rju (U), hence (i) is proved. □ 
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Lemma 8.1. The solution of 

(8.7) E«fo«)E ( l )(-iy(-j) k = B k u(x) 

1=0 j=0 V/ 

satisfies (3.2), and (3.5) respectively, where B := r) A. 

Proof. Suppose that the solution qi of (8.7) satisfies (I + l)qi+i = Bqi + 
Iqi, see (3.5), for < I < k, k > 0. Then (note that summations may be 
started at l,j = l for k > 0), 



fc+i l /A 

e <»(*.«) e [M-^i-jy 

1=1 7 = 1 W 



(8.8) C «b(x, u)E( •) ("l) J (-i) fc+1 = ^ fc+1 /«(x) 

transforms to 



k I / ,\ -, l-l 



EE (!• )(-l) i (-i) ,b+ 4n(S + r/)A,(x) = B*+ 1 A,( a 
z=ij=i V/ l >=o 

and after some straightforward algebra to 



k+i , k + r 



k+i /, 
%+1 (x,n)(-l) fc+1 E 

s fc+i /«(x) -ee ( l ) (-i) j (-j) k+l l ( b + r/ ) /«(*)■ 

1=1 j=l w 



Claim: For any fc > 
(8.9) 



(-i)*Ef*)(-i) i i* = *«- 

7=0 \-V 



.7 = 

This claim is easily proved by considering h(s) := (1 — e~ s ) k . On the one 
hand, 



— hi \ - — 

ds k [S) ~ ds k . 

3=0 



E(-l) 3 e- JS = E( fc )(-l) 3 (-j)^, 

7=0 7=0 V/ 
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and on the other, for \s\ small enough, h(s) = (s — ^s 2 + 



Therefore, 



ds k 



h(s) 



=v.=i:( k )(-iyH) k - 

=o j=0 W 



Using (8.9), (8.8) is equivalent to 

q k+1 (x, u)(k+ 1)1= B k+1 f u (x)- 



(_!)*+! ( l \-iyj k+ll M(B + rI)f u 

1 = 1 7 = 1 W l >=o 



(x) 



The rest follows from the next claim. 
Claim: For any k > and x £ 1 it holds 

k i /,\ 



(8.10) x 



k+1 



■D W EE 

1=1 j=l 



3 Ak+1 



l-l 



\{{X + 



r=0 



U(x + r). 



r=Q 



In order to prove this claim it is enough to show (8.10) for x = —k, —k + 
1, 0. Since for any natural to, 

1 Z-i 



n( 

r=0 



r — TO 



0. 



we have to show that for to < k, 



TO 



fc+1 



For = the above equality is obvious. Assume that we have proved the 
claim for k < n. Then it follows that 



I 

E 



-i) j j 



m—1 
TO ^ 
1=1 



TO 



s=0 



TO — 1 

, 1 , 
\ m—1 . 



)(-!)' 

TO — 1 



§G) ( - iyo+i) ' 



E 

i=i 



1) J J S 



n 



to (to — i) 

,=0 W 



S - TO™ +1 . 



□ 
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Proof of Theorem 3. 8. From the Taylor formula for semi-groups it follows 
that 



to kl rl Jo 

Due to (3.10), for < r < s < R u we have 

P T A r+l f u < sup ||P T || A r+l f u < sup ||P T || f-L + e y +1 (r + l)! 

for any e > 0. It thus follows that 



Psfu e ,,A k f u 



r+l 



s r+1 sup ||P T ||, 

0<r<s 



which converges to zero when r — ► oo, if \s\ < R u /(1 + sR u ). Since e > is 
arbitrary, the first statement is proved. 

The commutation property A k P t f u = P t A k f u and the boundedness of Pf 
for t > imply that for \s\ < R u , 



oo ,|fc 



OO life 



OO life 

S 



k=0 ' k=0 ' k=0 

Since P t f u G for all k > 0, (3.11) follows. 



< oo. 



Proof of Theorem 4.1. For r > define A r / U =: <7 r / u with f u {x) 
exp 



, and write 
A r+1 f u = A(g r exp(iu 1 

Leibniz formula implies 

A T+1 f u = E M*) E 

|q|>1 0<a 

E a «( x ) E 

|a|>l /3<a 

Hence, the following recurrent formula holds 
(8.11) g r+ i = E a <x( x ) E 



E a a (x)d xa ^ r exp(iit n 
H>i 



p\(a-[3)\ 



(5\{a-(3)\ 



d x f3g r d x c-p exp(iu T x) 
(\u) a ~Pd x (3g r exp(iu T x) 



H>1 /?<a 
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Similar formulas for derivatives of g T +i can be obtained: 

d xP g r+1 = V o U — mi {iu) a ~ 13 V -t-t-^ — rr<9^a Q d xP - v+0 g r . 

Since the underlying process is affine, all derivatives of a of order higher 
than one are zero and thus, by induction, g r is polynomial in x of degree at 
most equal r. We so get for \p\ < r + 1, 

d x pg r +i = ~1T~ — \T Yl mi ~ mi W^^Ca d xP - n +pg r 



By defining 



T r := max \d x /3g r \ , 
l/3|<r 



we obtain the following estimate for x E X, 

l<w +1 | < rva + iixid y. E E ^ra)i l " r "" D °»' 

v < p , 'l-yP 'V- \ a \>lP<a H ^ ' ' 
I»?I<1 

with |it| := [j tii | , |u n |]. Hence, by the simple relations 

\- 1 = (fc + ra- 1)! +fc 

fV, fc!(n-l)! " 

a,|a|=fc 

with ||ii|| = max^i^..^ we have 

(8.12) r r+1 < 2"r r (r + 2)(l + ||x||)^2 fc (l + ||n||)^ 

k>l 

= 2»r r (r + 2)(l + ||x||)0(||u||), 

where the series in (8.12) is convergent due to assumption (4.2). As a con- 
sequence, (4.3) holds. □ 
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Proof of Proposition 4-6- From (8.11) we have with do := 0, 
9r+i = E^ a <* i3\(a- p)\ W*~ Pgr " ( 7 1 py ^^r 1 ^* 1 ^ 

= E 9r,l+f3\ + n^ )^ 7 l| 7 +/3|<r^/3 
/3, 7 >0 \ P / 

= E- 7 E ^(VV? 

|7|<r \P\<r-\l\ \ 1 / 

+ E E ^f 7 /) E *k^ +K 

\y\<r |/9|<r-|7| \ M / /t, |/e|=l 

= E - 7 E ^(Vk 

|7|<r+l |/9|<r-|7| \ ^ / 

+ E x7 E E ^, 7 - K +/3( 7 ^ +/ j^, K , 



H<r+1 |«|=1, K<7/3, |j3|<r+l-|7| 



where empty sums are to be interpret as zero. The second recursion follows 
from (r+l)h r +i = r) h r +i+rh r with A(h r f u ) = h r+ if u and h r+ i computed 
via (4.6) with g r replaced by h r . □ 

References. 

[1] Amadori, A. L. (2003) Nonilinear integro-differential evolution problems arising in 
option pricing: a viscosity solutions approach, Differential and Integral Equations 16, 
787-811. 

[2] Dieudonne, J. (1960). Foundations of Modern Analysis. Academic Press, New York. 
[3] 0KSENDAL, B. and Sulem, A. (2007). Applied Stochastic Control of Jump Diffusions. 
Springer. 

[4] Duffie, D., Pan, J. and Singleton, K. (2000). Transform analysis and asset pricing 

for affine jump diffusions. Econometrica, 68, 1343-1376. 
[5] Duffie, D., Filipovic, D. and Schachermayer, W. (2003). Affine processes and 

applications in finance. Annals of Applied Prob., 13, 984-1053. 
[6] Filipovic, D. (2001). A general chracterization of one-factor affine term structure 

models. Finance Stoch., 5,389-412. 
[7] Eberlein, E. and Papapantoleon, A. (2005). Symmetries and pricing of exotic op- 
tions in Levy models. In Exotic option pricing and advanced Levy models, A. Kypri- 

anou, W. Schoutens, P. Wilmott (Eds.), Wiley. 
[8] Glasserman, P. and Kyoung-kuk Kim (2007). Moment Explosions and Stationary 

Distributions in Affine Diffusion Models. Working Paper. Columbia Business School. 
[9] Kampen, J. (2006). The WKB-Expansion of the fundamental solution 

of linear parabolic equations and its applications, book manuscript. At 

http://papers.ssrn. com/abstract=914122, Soc. Sc. Res. Net. 



imsart-aap ver. 2007/12/10 file: bks_annal_prob.tex date: July 8, 2008 



30 



D. BELOMESTNY ET AL. 



[10] Keller-Ressel, M. and Steiner, Th. (2008). Yield Curve Shapes and the Asymp- 
totic Short Rate Distribution in Amne One-Factor Models. To appear in Finance 
Stock. 

[11] Protter, P. (1990). Stochastic Integration and Differential Equations. Springer. 
[12] Sato, K. (1999). Levy Processes and Infinitely Divisible Distributions. Cambridge 
University Press. 

Weierstrass Institute for Applied Analysis and Stochastics 

mohrenstr. 39 

10117, Berlin, Germany 

E-mail: belomest@wias-berlin.de 

kampen@wias-berlin.de 

schoenma@wias-berlin.de 



imsart-aap ver. 2007/12/10 file: bks_annal_prob.tex date: July 8, 2008 



